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Summary 

A simple energy approach to study the problem of control -structure inter- 
actions in large space-structures is presented. For the illustrative case of a 
free-free beam, the vibrational energy imparted during operation of constant, 
step, and pulsed thrusters is found in a non-dimensional closed form. Then 
based on a parametric study, suggestions are made on the choice of parameters 
to minimize the control -structure interactions. The study of this simple sys- 
tem provides physical insight and understpnding for more complex systems. 
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Notation 


E.,E energy and its non-dimensional form 

' flexure rigidity coefficient 

force due to thruster 
transverse moment of inertia 
Ri ratio of energy in flexure, flexure mode I 

kinetic energy 
strain energy 
mode-shape 
rZi/Ti 
pulse number 

roots of eq. (9) defining mode-shape 
length > 

mass per unit length 
total number of pulses 
frequency 
time 

distance from one end of beam 
transverse deflection 


oo 



time between pulses 
period of t mode 
pulse-width 
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Subscript 


f 

I 

max 

m’m 

r 


flexure mode 
mode 


maximum 
mi n imum 
rigid mode 

Dots indicate differentiation with respect to time and primes indicate differen- 
tial with respect to X , 


Introduction 

The problem of control -structure interactions in large spacecraft has been 
of concern in recent years. A number of studies have been conducted on the sub- 
ject by many authors. Generally one has to resort to extensive numerical simu- 
lation. Any a priori knowledge of the system helps in reducing the computational 
effort involved. This paper presents a simple energy approach which provides 
considerable insight into the problem. As an illustration a free-free uniform 
beam subjected to continuous, step and pulsed operation of a transverse thruster 
is chosen. It represents a rectangular solar power satellite or a slender 
rocket with large length to width/diameter ratio. To control the orientation 
(rigid-mode), transverse thrusters, usually placed near the ends, are operated. 
Their operation can impart significant undesirable vibrations. 

Here the ratio of energy in flexure mode and rigid mode is found in a closed 
form. It is then used for an extensive parametric study. Non-dimensionalization 
makes the results applicable to a wide range of systems. Based on the analysis, 
suggestions are made for a strategy minimizing the control -structure interac- 
tions. 
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Analysis 

Consider a free-free beam (Fig. la) subjected to a concentrated transverse 
load (thrust) F. Under the action of the force the beam will translate, rotate 
and vibrate. The first two are designated as rigid-mode, and the last is called 
flexure-mode of motion. Let us consider three types of loading - constant, step 
and pulsed. 


A. Constant (continuous) Thrust 

Ion-thrusters generate a loading of this nature (Fig. lb). Under its ac- 
tion, energy in the rigid mode is given by 


Er - Tr =(I co^+ m{ V^)/z 


( 1 ) 


where CO and D are rational and translational speed. If the thruster is at one 
end 0 we get 


Er = XF^byrnl 


( 2 ) 


The kinetic and strain energy in the flexure mode are given by 

If = 


- -i I Ei(yfdx 


(3) 


Under concentrated loading the response of a beam can be written in terms of 
modal components as^ 


M = I K (t -t') d b' 

I i ^ o 


( 4 ) 
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whereX({ is mode calculated at the point of application of F , and^^*^. 
For a constant force, the response is 

(5) 

' t 

For a free-free beam the normal modes can be taken as 

X= ± (coskx4cosli k*')±c(sinka; + sinlikx) 


where . 

C - |(sinh k5 4- •Sinkf)/(sink k5 — SinkI) j ^ 


(i. =>/il/nnl<L^ 

coskil cosh kl = i 


(7) 

( 8 ) 
(9) 


Noting that 

■jV'fdo: = 5 


; J dx = 0 i=^j 


(10; 


and using (5) - (8), eq. (3) can be written as 

T y V?. - sin^ Pc^ 


p/ 


u 


f 3mi 


7 . Xfi (i- cosp.iryp.2 


( 11 ) 
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The total energy in the flexure mode is therefore 

g, . Z ( i _ cos Kt)/Pi 


■f 


(12) 


For the force acting at one end X^£=:4 » and 

Ef = PibA)/p- 


(13) 


Eq. (2) and (13) can be written to yield non-dimensionalized components of 
energy as 


E/ = Er mi/2F"xf = 

EJ = Efnal/^F^li" = l(Tt/rTi)’"^in>t/xO 


(14) 


(15) 


The energy-ratio is given by 

K= Ef/Er = 1 (-S'*! / ^0^ = 2^1 


(16) 


where “ 7T ^/~C^ 

Thus we obtain simple expressions for flexure-energy and energy-«*atio. A high 
value of R corresponds to large control -structure interaction. Our aim is to 
minimize it. This will be discussed in the next section. 


B. Step Thrust: 

Generally! chemical thrusters are operated for a limited time like a step- 
function (Fig. Ic). During the application of force all the relations found 
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above hold. After the termination of thrust at . the energy will remain con 
stant and will correspond to the value at Hence witht replaced by fcj_ , 
relations (12) - (16) are applicable. The response after tj^is given by 

1'=°^ cos|=.t}/p> ( 17 ) 


C. Pulsed Thrust: 

In many applications control is achieved using pulsed thrusters. In this 
case control -structure interactions can severe if the system parameters, namely 
thrust level F , pulse-width S , period between pulses "C , and number of pulses 
n , are not chosen properly. 

The response at the end of n uniform- pulses can be expressed as 


y = I Pc tjc)] (la, 

where tjf and tjo refer to the final and initial time of pulse. They can 
be expressed as 


t3f= (j-i)I+5 


tjo = 0-i)T 


(19) 


Using (19), the J series in (18) can be summed to yield 

a mu ^ Sin PiT/i L ^ ■' 


( 20 ) 
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Using (20) in (3) with (6) and (10) yields the vibrational kinetic and strain 
energy at the end of n p ses as, 

Tf z unYcC^-^ix)- s',n PcCb- 

J ^nu fSj^Sin P^T/2. i 


U 


f 


« Z "f i) 


( 21 ) 


The total flexure-mode energy Is 


CospiCt-n^T)}] 


/Xo£ln^ 

ttiK ^ V Pi -Sifl- ^ ^ 


(22) 


With the thruster at one end, this can be written in non-dimensional form as 


C* ^ 


J Xi.-Sin nrrryxc -Sin iicS^/Xi 7^ 

TTT^ (ft/Xc J 


(23) 


Corresponding component of energy in the rigid-mode is given by 


- ("Vr,)’ 


(24) 


From these, "R can also be written. Note the similarity with (16) - (18). The 
limiting values of the energy introduced into the flexure-mode can be easily 
found from (23) as follows. 


Since 



oir x/tc = 

n x/Xc • * and x/lt ¥ ‘hi-e^er 

(25) 
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we have from (23), 




i ‘min ** ^ 


(26) 


Also, 

;f 


=. E 


C - s Q 

jc mav. 


at* • ^/Xi - (N 4- i/z.') TT 

QS -*“CO 

a(- cT/tc = Vl 


(27) 


where H is an integer. With a proper choice of ,T1 it is possible to 

minimize the total , as discussed in the next section. 


Parametric Study 

In the case of continuous constant thrust, the only parameter under control 
is the thrust level F . Both Ef and Ep increase as square of F . The energy- 
ratio‘s varies with time. Initially for each mode, is near unity. It goes 
down to zero when "t equals the period of that mode (Fig. 2a). Fig. 2b compares 
the energy-ratio in the first 10 modes at three instants from the start. Gen- 
erally the first three modes dominate. At some instants the higher modes may 
contain more energy than the lower ones’. The overall energy ratio is very high 
in the beginning. Gradually, it reduces with the increase in time (Fig. 2c). 
Fig. 2d shows that initially only flexure modes are excited. Then the rigid- 
mode gains energy and keeps building up. Energy in the flexure-mode, hov/ever, 
keeps fluctuating between 0.01 and 0.11, which occur near NF" and 7T , 

respectively (Fig. 2d). 
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These observations are applicable for the case of step-thrust also. Here, 
for a given amount of control ( Ef constant), it is possible to trade between F 

and to minimize the energy introduced into the flexure^mode using Fig. 2d. 

Generally a choice of t^^nTTTj^ will give minimum flexure-energy. After 
choosing , F can be fixed. From eq. (13) and (15) one may also note that 
for a given and F more modes are excited for a more flexible system. 

Use of pulsed thrust introduces 4 parameters, F , cT , "C , and n , any 

three of which can be chosen at will for a given amount of control. Fig. 3 and 

4 generated using (23) show the effect of varying these parameters. Fig. 3(a) 
shov/s the maximum energy in a mode (Eq. 26) as a function of the pulse-width. 

As expected, for a pulse-width equal to an intepral multiple of the modal period, 
the energy in that mode is minimum. However, this may enhance energy in the 

other modes. So the choice of the pulse-width which minimizes the total flexure 

energy will require some effort and t-ade-offs. For a given thrust level and 
rigid-mode energy. Fig. 3(b) shows that the maximum energy in a mode increases 
as we increase the number of pulses (and reduce the pulse-width). They asymp- 
totically reach the value corresponding to the rigid-mode. The energy in higher 
modes fluctuates between zero and a small value up to some n before increasing 
monotonically. Hence it seems ^ood to keep n low. It is possible that energy 
in some modes may be more than in the lower modes. It is interesting to note 

that for the case considered only the first mode is excited at n= 3. 

After fixing and n (and thereby, F for a given E^) to yield the mini- 
mum of peak flexure-energy, it is possible to then reduce it further by a proper 
choice of the pulse-period Z . Generally it is desirable to keep Z as small 
as possible to complete the operation soon. However, one should avoid resonance 
with any mode, otherwise, the energy in flexure-mode can be very high. This is 
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shown on Fig. 4, which presents variation of ^ as a function of X. for two 
cases ( n- 4 ^ <5 Xx and o-i ) with a given level of 

thrust and rigid-mode e- argy. With the increase in the number of pulses higher 
modes are also excited and the resonance peaks get higher and narrower. The 
choice of T . therefore, becomes more critical. At (’C/Xi )^*»^the pulsed 
thrust case reduces to a step function. Fig. 4(c) shows the total energy in 
flexure-mode as a function of . For n *4 , x/XjsO-^ results in min- 

imum total energy. For Hr-io minima occur atx/tx " 

At higher r\ a slight change from the 'optimum X* can result in a significant 
increase in . 

Based on these observations, a possible approach to minimize the control- 
structure interactions can be as follows: 

• from attitude -control considerations, select the mode of control. 
This will also fix the energy per operation in the rigid-mode. 

• thrust level is decided by the available hardware 

• for the given Ej- and F generate plots like Fig. 3 and select 6* 

and n to give lowest p; * 

r mQ-)^ 

• for the given Ej- , F and n use eq. (23) to obtain a plot like 

Fig. 4(c). Then choose to yield the lowest 

Note that one may interchange the last two steps. 

Alternatively, one may employ the variation of parameter approach to find 

the set of parameters corresponding to minimum r* (eq. 23). For this purpose 

f 

the series may be truncated judiciously. 

After establishing , the maximum amplitude of vibration in a mode can 
be found by equating to the component of strain-energy in that mode. 
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Co ncluding Remarks 

A simple approach to study the problem of control -structure intera:t1w . 
in large space-structure is presented. The energy considerations applied to 
the case of a free-free beam, subjected to continuous, step and pulsed thrust, 
yield much insight into the system behavior. This also leads to find a way to 
minimize the interactions. The approach can be used for many systems. It should 
reduce the computational effort involved in such analysis significantly. 
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(d) Pulsed thrust 


FIGURE 1: 


A FREE BEAM SUBJECTED TO CONCENTRATED LOAD 










